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Abstract The dynamics of gradient and Hamiltonian flows with particular application to flows on adjoint orbits of a 
Lie group and the extension of this setting to flows on a loop group are discussed. Different types of gradient flows that 
arise from different metrics including the so-called normal metric on adjoint orbits of a Lie group and the Kahler metric 
are compared. It is discussed how a Kahler metric can arise from a complex structure induced by the Hilbert transform. 
Hybrid and metriplectic flows that arise when one has both Hamiltonian and gradient components are examined. A 
class of metriplectic systems that is generated by completely antisymmetric triple brackets is described and for finite- 
dimensional systems given a Lie algebraic interpretation. A variety of explicit examples of the several types of flows 
are given. 
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1 Introduction 



Dynamical systems, finite or infinite, that describe physical phenomena typically have parts that are in some sense 
Hamiltonian and parts that can be recognized as dissipative, with the Hamiltonian part being generated by a Pois- 
son bracket and the dissipative part being some kind of gradient flow. The description of Hamiltonian systems has 
received much attention over nearly two centuries and, although some forms of dissipation have received general at- 
tention, the understanding and classification of dissipative dynamics is a much broader topic a nd consequently less 
well developed. Early moder n treatments of geometric Hamiltonian mechanics include those of ISouriaul 1 197011 and 
Abraham and MarsdenI 1 1978 1 , and the literature on this top ic is now immen se. A spe cial type of gr a dient flow tha t pre- 
serves invariants, the double bracket formalism d escribed in BrockettI 1 1991 ] (see, e.g., Bloch 1 1990l| . Bloch 1 2003 1 ), is a 

r 1! ^1- _i !„ ! _i r i i_ / Ini i- T7-- _i j n T j j r» i! i r\r\ yi' ' i r\r\/rh\ i : ii 



formalism that occurs in a variety of contexts (see'Bloch, Krishnaprasad, Marsden, and Ratiu f 1994, 1996']) and is well- 
adapted to practical numerical computations (see iVallis. Carnevale and Young. [1989. 1: ,Flierl and Morrison 1.201 111 ). Ex- 
amples of infinit e-dimensional gr a dient flows inc lude the Cahn-Hilliard systems (see lOttd 11200 ill ) and the celebrated 
Ricci flows (see HamiltonI 1 1982 1: Chow 1 2004 1). which are nonlinear dif fusion-like equations. A general form for 
combined Hamilto nian and g r adient flows was describ ed in lMorrisonlil986ll. where such flows were termed metriplec- 
tic flows (see also Oettingei 1 2006ll : Morrison 1 2009ll : Liero and Mielke 12012 1). Thus, it is evident that there are a 
variety Hamiltonian and dissipative flows, and the purpose of this paper is to explore the form and geometric structure 
of such flows in both the ode and pde contexts. 
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Specifically, in this paper we discuss the dynamics of gradient and Hamiltonian flows, with particular application to 
flows on adjoint orbits of a Lie group and the extension of this setting to flows on a loop group. We compare the different 
types of gradient flows that arise from different metrics, in particular, the so-called normal metric on adjoint orbits of 
a Lie group and the Kahler metric. We discuss how a Kahler metric can arise from the complex structure induced 
from the Hilbert transform. We also consider flows that arise when one has both Hamiltonian and gradient structures 
present. In particular, we discuss metriplectic flows, flows that produce entropy while conserving energy. We consider 
such flows in both the finite and infinite settings, and discuss a general class of metriplectic flows that arise from 
completely antisymmetric triple brackets. For finite systems, we show how the triple bracket has a natural Lie algebraic 
formulation, and for infinite systems we give a procedure for constructing a quite general class of metriplectic pdes. 
We also consider, hybrid flows, of Hamiltonian and gradient form, that dissipate energy. Several examples of hybrid 
and metriplectic flows are given, including finite systems such as the Toda lattice on K. and metriplectic so(3) brackets. 
Various infinite-dimensional example s including a 1 + 1 dis sipative systems that conserves energy, and hybrid systems 
such as the KdV with dissipation, the lOtt and Sudani il969ll equation that describes Landau damping, and others. 

The paper is organized as follows. In section |2] we review material need for latter development. In particular, we 
discuss metrics on adjoint orbits, Toda flows and the double bracket formulation. Sections[3]and|4]contain the main new 
results of the paper as described above. In section|3]we discuss metrics on loop groups and related gradient flows, while 
in section|4]we discuss our results on metriplectic systems, in both finite- and infinite-dimensions, and give examples. 



2 Metrics on adjoint orbits of compact Lie groups and associated dynamical systems 
2.1 Double bracket systems 

Let Qi, be the compact real form of a complex semisimple Lie algebra g, G„ a compact connected real Lie group with 
Lie algebra g„, and K the Killing form (on g or g„, depending on the conte xt). 

The "normal" metric on the adjoint orbit ff of Gj, through Lq G Qu (see Atiyahl 1 1982 1. I Besse . 2008 . Chapter 8]) is 
given as follows. Decompose orthogonally g„ = 0^ ® Qul, relative to to the invariant inner product ( , ):=—«■( , ), 
where QuL '■= kerad/, is the centralizer of L and = rangead/,; as usual, adi := [L, ■]. For X G g„ denote by £ g^ 
and Xl S Qul the orthogonal projections of X on g^ and Qul, respectively. Recall that a general vector tangent at L 
to the adjoint orbit ^ is necessarily of the form [L,X] for some X G 0„. The normal metric on ^ is the G[,-invariant 
Riemannian metric given by 



normal 



-{X^Y^-) 



forany X,y G g„. 

Fix N G Qu and consider the flow on the adjoint orbit ^ of G„ through Lq G Qu given by 



-Lit):=[Lit),[Lit),N]], 



L(0) = Lo G g„ . 



(1) 



(2) 



We recall the following well-known result ( Brockett J 199 1 1 . [Brocket 3 lll994ll. iBloch, Brockett. and Ratiul il990ll . 
Bloch. Brockett. and Ratiul il992ll . lB"loch. Flaschka. and Ratiul j 1 990ll . lBloch and IserlesI 1200511 ). 



Proposition 1. The vector field given by the ordinary differential equation (|2]l is the gradient of the fimction H{L) = 
k{L,N) relative to the normal metric on 6. 

Proof. By the definition of the gradient grad//(L) G Tiff C g„ relative to the normal metric, we have for any L G ^ 
and 5L g g„, 

dH{L) ■ [L,8L] = (grad//(L), [L, 5L])nonnai (3) 

where • denotes the natural pairing between 1 -forms and tangent vectors and [L, 5L\ is an arbitrary tangent vector at L 
to ff. Set grad//(L) = [L,X] = [L,X^]. Then OJ becomes 

-([L,5L],A?) = ([L,X],[L,5L])„ormal 



or, equivalently, 



Gradient flows and metriplectic systems 3 

([L,A^],5L) = (X^5L^) = (X^5L). 
Since [L,N] e fl^, this implies that = [L,N], and hence grad H{L) = [L, [L,N]], as stated. ■ 

The same computation, for a general function H G C°°(gi,), yields 

gradi/(L) = -[L,[L,V//(L)]] (4) 

where VH{L) denotes the gradient of the function H relative to the invariant inner product ( , ) := — K"( , ), i.e., 
dH{L) X = (V//(L),X) for any X e 0„. 



2.2 The finite Toda system 



The double bracket equation (|2]i is intimately related to the finite non-compact Toda lattice system. This is a Hamilto- 
nian system modeling n particles moving freely on the x-axis and interacting under an exponential potential. Denoting 
the position of the A:th particle by xj^, the Hamiltonian is given by 

^/t=l k=l 



and hence the associated Hamiltonian equations are 

dH dH 

oyk oxk 



(5) 



where we use the conventions e^o^^i = g^n-^n+i = Q, which corresponds to formally setting xq = —°o and x„+i = +0°. 

This s ystem of equatio ns has an extraordinarily rich structure. Part of this is revealed by Flaschka's change of 
variables ( Flaschkal 1 1974 ]) given by 



fli=^e^'^*"'^*+''/^ and fei = -^y^. 

which transform ^ to 

{ ak=ak{bk+{-bk), ^=l,...,n-l, 
I^A: =2(flf-af_i), k=\,...,n, 
with the boundary conditions aq = a„ = 0. This system is equivalent to the Lax equation 



(6) 



—L^\B,L]=BL-LB, 
dt 



where 



fll •■■ \ 
a\ ■ ■ ■ 



L = 



■■■ fl„_i 
\0 ••• fl„_i b„ J 

If 0{t) is the orthogonal matrix solving the equation 



B = 



/ fli 

— fli fl2 



V ••• 



\ 




fl„_i 
-fl„_i / 



(7) 



(8) 



dt 



= B0, 0(0) = Identity, 



then from (|7]i we have 
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Thus, O LO =L(0), i.e., L{t) is related to L(0) by conjugation with an orthogonal matrix and thus the eigenvalues 
of L, which are real and distinct, are preserved along the flow. This is enough to show that this system is explicitly 

solvable or integrable. Equivalently, after fixing the center of mass, i.e., setting /?! H ^b„~Q, the « — 1 integrals in 

involution whose differentials are linearly independent on an open dense set of phase space { (a i a„_ i , i ,...,/?„ ) 
^1 + • • • +^« = 0} are TrL^, . . . ,TrL". 



2.3 Lie algebra integrability of the Toda system 

Let us quickly recall the well-known Lie algebraic approach to integrability of the Toda lattice. Let g be a Lie algebra 
with an invariant non-degenerate bilinear symmetric form (,), i.e., {[^,r\],Q = (<^j['?jC]) for all £,,ri,^ e g and 
{£,,■) ~0 implies £, = 0- Suppose that t,s C g are Lie subalgebras and that, as vector spaces, g = 6©s. Let TTf : g — > 6, 
TTs ; g s be the two projections induced by this vector space direct sum decomposition. Since g 9 i ') G 0* 

is a vector space isomorphism, it naturally induces the isomorphisms = s*, ^ ^ * gy non-degeneracy of ( , ), we 
have g = 5^ © t^; denote by k^l : g — ?► t^, n^± : g — >■ s-'- the two projections induced by this vector space direct sum 
decomposition. In particular, g, s^, all carry natural Lie-Poisson structures. The (-)Lie-Poisson bracket of s* = 
is given by 

{cp,wm = -{^,[n,Vcp{^),n,Vvm, ^et^, (9) 

where ^, \/a : — > R are any smooth functions, extended arbitrarily to smooth functions, also denoted by (p and on 
g and V^, Vi// are the gradients of these arbitrary extensions relative to ( , ). This formula follows from the fact that the 
gradient on of ^|{±, which is an element of s due to the isomorphism = s*, equals TisVcp. Thus, the Hamiltonian 
vector field of i// e C°°({^), given hy <p = {(p , Xj/} for any (p e C°°(t^), has the expression 

X^(^) = -7rjx[7r.Vv/(^),^], <^et^ (10) 

with the same conventions as above. 

If G C°°(g) is invariant, i.e., [Vv/(Q, Q—0 for all C S g, then ( fTot simpUfies to 

X^{^)^[7ZtVw{^U] = -[7^.^V{^U], ^Gfi^- (11) 



The Adler-Kostant-Symes Theorem (see lAdleJ il979ll . lKostand il979ll .ls vmes I lll980allbll . and|Ratiu| il980ll for many 



theorems of the same type) states that if (p and y/ are both invariant functions on g, then {(p, = on 6^ which is 
equivalent to the commutation of the flows of the Hamiltonian vector fields (fTTT i. 

Suppose that G = KS, where G is a Lie group with Lie algebra g and A',5' C G are closed subgroups with Lie 
algebras i and s, respectively. The writing G = KS means that each element g G G can be uniquely decomposed as 
g = ks, where k G K and s € S and that this decomposition defines a smooth diffeomorphism K xS ^G. The coadjoint 
action of S on s* has the following expression, if 0* is identified with 1^ via {,): if s <E S, ^ G 4^, then s- £, ^ n^i_ Ad, ^ , 
where Adj ^ is the adjoint action in G of the element ^ G 5 C G on G C g. 

For the Toda lattice this general setup applies in the following way. Let G = GL(n,M), K = SO(n), S = 
{invertible lower triangular matrices}, G = KS is the Gram-Schmidt orthonormalization process, g = g[(«,M), t = 
so(n), s = {lower triangular matrices}, , rj) := Tr(^77) for all ^ , T] G g[(n,M), 6^ = sX)m{n) the vector space of sym- 
metric matrices, and = n, the nilpotent Lie algebra of strictly lower triangular matrices. The set of matrices L in 
([8]) is a union of S'-coadjoint orbits parametrized by the value of the trace; for example, the set of trace zero matrices 
L of the form ^ equals the 5-coadjoint orbit through the symmetric matrix that has everywhere zero entries with the 
exception of the upper and lower first diagonals where all entries are equal to one. Thus, the Toda lattice is a Poisson 
system whose restriction to a symplectic leaf is a classical Hamiltonian system with « — 1 degrees of freedom. The 
Hamiltonian of the Toda lattice is 5 TrL^ and the fk{L) := j TtL'^, k= 1 ,...,« — 1 are the « — 1 integrals in involution 
(by the Adler-Kostant-Symes Theorem) and are generically independent. 
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2.4 The Toda system as a double bracket equation 



If is the matrix d iag|l,2, . . . ,» |, the Toda equations (|7J may be written in the double bracket form (|2]i forB :— [N,L\. 
T his was shown in Bloch J 1990ll: the c onsequences of this fact we r e furth er an alyzed for general compact Lie alge bras 
in iBloch. Brocket!, and Ratiul il990ll . iBloch. Brockett. and Ratiul lll992ll . and iBloch. Flaschka. and Ratiul il990ll . As 



shown in Proposition [T] the double bracket equation, with L replaced by iL and by iA^, restricted to a level set of the 
integrals described above, i.e., restric ted to a generic adjoint orbit of SU (n), is the gradient flow of the function TrLA^ 
with respect to the normal metric; see lBloch. Flaschka. and Ratiu il990ll for this approach. 

This observation easily implies that the flow tends asymptotically to a diago nal matrix wit h the eigenvalu es of 
L(0) on the diagonal and ordered according to magnitude, recovering the result of lMosen 111 97511 . ISvmesI lll982ll . and 
Deift. Nanda. and Tomeillll983ll . 



2.5 Riemannian metrics on 



Now, we recall that, in addition to the normal metric on an adjoint or bit, there are o ther natural G„-invarian t metrics: 
the induced and the group invariant Kahler metrics (as discussed in I AtiyahL 1982 , §4], Atiyah and Pressley. I,1983il . 
and iBesseLl2008l Chapter 8]). 

Firstly, there is the induced metric b on ff, defined by b := i* {~k{ , )), where i : ^ g„ is the inclusion and 
( , ) := — k:( , ) is thought of as a constant Riemannian metric on g„. Therefore, 



b{mL,XULJ]) -.^ {[L,XUL,Y]) 



(12) 



for any L € if,X,Y(i g,,. The induced metric on is also G„-invariant. 

Secondly, there are the G„-invariant Kahler metrics on & compatible with the natural complex structure (of course, 
induced by the complex structure of G). These are in bijective correspondence (by the transgression homomorphism) 
with the set of G„-invariant sections of the trivial vector bu ndle over w hose fiber at L G ^ is the center of ker (ad/^) 
and whose scalar product with all positive roots is positive ( lBessell2008l Proposition 8.83]). Among these, there is the 
Ga-invariant Kahler metric b2 which is compatible with both the natural complex structure on and has as imaginary 
part the orbit symplectic structure; ^2 is called the standard Kahler metric on ff. 

The G„-invariant Riemannian metrics on a maximal dimensional orbit & are completely determined by T-invariant 
inner products on the direct sum of the two dimensional root spaces of g„, which is the tangent space to at the point 
Lo G t in the interior of the positive Weyl chamber; recall that intersects the positive Weyl chamber in a unique point. 
The negative of the Killing form induces on each such 2-dimensional space an inner product. This inner product, left 
translated at all points of by elements of G„, yields the normal metric on ff. Any other G„-invariant inner product 
on & is obtained by left translating at all points of the inner product on this direct sum of 2-dimensional root spaces 
obtained by multiplying in each 2-dimensional summand the inner product with a positive real constant. 

Since Lq lies in the interior of the positive Weyl chamber (because is maximal dimensional), o:(Lo) > for all 
positive roots a of Qu- Then the constant by which the natural inner product on the 2-dimensional root space needs 
t o be muhiplie d in order to get the standard Kahler metric is a(Lo), w hereas to get the induced m e tric, i t is Oc(Lo)^ 
( llAtivah[[l982i Remark 2 in §4]). We can formulate this differently, as in lBloch. Flaschka. and Ratiul Ill990ll . Since, by 
dUli and O, 

b{L){mxULJ]) = {[L,XULJ]) = {[L,X\[LJ^]) = {-[L, mx%Y^) = {-[L, [L,X^]]\y^) 
= {-^Al%Xl%Y])_^ 

b{L){[L,XUL,Y]) = biiL){£/{L)^[L,X],[L,Y]), (13) 



we have 



where we denote now by bi the normal metric and jz/ (L) := ij (iad^)^ is the positive square root of (iad/,)^ = — ad|^ = 
£/{L)'^. The standard Kahler metric on ^ is then given by 



b2{L)[L,X], [L,F]) = bi{s^{L)mx],[L,Y]). 



(14) 
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Note that, as opposed to the normal and induced metrics which have explicit expressions, the standard Kahler metric 
on requires the spectral decomposition of (L) at any point 6 . Or, as explained above, one expresses it at the 
point Lo in the positive Weyl chamber in terms of the positive roots and then left translates the resulting inner product 
at any point of & . The normal metric does not depend on the operators (L), whereas the standard Kahler and induced 
metrics do. 



3 Gradient flows on the loop group of the circle 

In this section we introduce three weak Riemannian metrics on the subgroup of average zero functions of the connected 
component of the loop group L(5'') of the circle, analogous to the normal, standard Kahler, and induced metrics on 
adjoint orbits of compact semisimple Lie groups. Of course, we shall not work on adjoint orbits of this group because 
they degenerate to points, L(5'') being a commutative group. Then we shall compute the gradient flows for these three 
metrics. 



3.1 The loop group ofS^ 



Recall (e.g., Presslev and Segal il986ll ) that the loop group L(5'') of the circle 5' consists of smooth maps of 5' to 
S^. With pointwise multiplication, L(5'') is a commutative group. Often, elements of L(5'') are written as e^^, where 
/ G L(]R) := {g : [—7t, tt] — > M | g is C°°, g{n) = Tt) + 2nn, for some « G Z}; « is the winding number of the closed 
curve [— ;r, tt] 9 f i— ?• e^^^'^ G 5' about the origin. More precisely, there is an exact sequence of groups 







L(I 
Inn; 



exp 



/ 



L(5') 



Z — > 

f(n)-f(-iz) 



which shows that kerexp = Z and cokerexp = {0}. Thus the connected components of L(5') are indexed by the 
winding number The connected component of the identity L(5'' )o consists of loops with winding number zero about 
the origin. 

If one insists on working with smooth loops, then o ne can consider L(5'') and L(5Mo as Frech et Lie groups either 
in the convenient calculus of lKriegl and Michorl 1199711 or in the tame category of lHamiltonI 1198211 . 

Alternatively, one can work with loops e'^ for / : \—n, 7r ] — > R of Sobolev class H\ where ^ > 1 (or appropriate 



W^'P or Holder spaces). By standard th eory (see, e .g..lPalaia Ill968 | or [ Adams and Fournierl rt2003ll X it is checked that 
L(5') is a Hilbert Lie group (see, e.g.. lBourbaM ll97lll or lNeebl l20()4l '). We shall not add the index s on L(R) and 
L(5'); from now on we work exclusively in this category of W Sobolev class maps and loops. A simple proof of the 
fact that L(M) is a Hilbert Lie group was given to us by K.-H. Neeb. First, note that L(R) is a closed additive subgroup 
of the Hilbert space H'{M) {/i : M M | /i of class H'}. Second, L(R) = L(R)o x Z as topological groups, where 
L(M)o {g G L(]R) | g{n) — g{~n)} is the closed vector subspace of //'(M) consisting of periodic functions; hence 
it is an additive Hilbert Lie group. Therefore, there is a unique Hilbert Lie group structure on L(M) for which L(R)o 
is th e connected component of the identity. For general criteria that characterize Lie subgroups in infinite dimensions, 
see llNeebll2006l Theorem IV.3.3] (even for certain classes of Lie groups modeled on locally convex spaces). Third, 



since exp : L(IR) L(5'^ ) maps bijectively each connected component of L(R) to a connected component of L(5'^ ), it 
induces a Hilbert Lie group structure on L(5'). 

The commutative Hilbert Lie algebra of L(5'') is clearly //*'(5',M) := {m : 5' — M | m of class the space of peri- 
odic//' maps, and the exponential map exp : //•*'(5'\M) — L(5'') is given by exp(M)(0) = e"'(^), where Q G R/27rZ = S^. 



3.2 The based loop group ofS^ 

The inner product on the Hilbert space (5' ) of real valued functions on 5' is defined by 
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2kJ- 



Following PresslevI il982ll and lAtivah and Presslev 1 1983 1. we introduce the closed Hilbert Lie subgroup L(5') : = 
{(j? eL(5') I (p(l) = l}ofL(5') whose closed commutative Hilbert Lie algebra is L(M) := {ueH'{S\R) \ m(1) =0}. 
The exponential map exp : L(R) 3 e"' G L(5^ ) is a Lie group isomorphism (with L( K) thought of as a commu tative 
group relative to addition), a fact that will play a very important role later on (see also jPresslev and Segal , Il986l page 
151, §8.9]). 

There is a natural 2-cocycle co on L(M), namely 



1 



d0M'(0)v(0) = (m',v) 



2n J-Ti 

where u' := du/dQ. Therefore, there is a central extension of Lie algebras 

— > R — > L(R) — > L(M) — > 



(15) 



which, as shown in lSegall 11198111 . integrates to a central extension of Lie groups 



1 



■5' 



■L(5i 



L(5i 



L 



The "geometric duals" of L(K) and L(M) — R ® L(R) are themselves, relative to the weak L^-pairing. It turns out that 

the coadjoint action of L(5'') on L(R) preserves {1} ©L(R) so that, as usual, the coadjoint action of L(5'') on L(R) is 
an affine action which, in this case, because the group is commutative, equals 



AdVM 



L 
f 



(log|/|)' e'feLis'), neLi 



Thus, the orbit of the constant function is L(5'')/5'' (where the denominator is thought of as constant loops), i.e., it 
equals L(5''). Therefore, every element u E L(R) of its Lie algebra has, in Fourier representation, vanishing zero order 
Fourier coefficient , i.e., m(0) — 0. 

Thus, the based loop group is a coadjoint orbit of its natural central extension and, a c cording to jlH has three dis - 
tinguished weak Riemannian metrics. These were computed expUcitly in lPressleyl lll982ll . lAtiyah and PresslevI 11 198311 . 



Pressley and Segal 1.1 9861 : we recall them below. 



3.3 L(5^) as a weak Kdhler manifold 

Note that on L(R), the cocycle ( fTSl ) is weakly non-degenerate. Therefore, left (or right) translating it at every point of 
the group L(5') yields a weakly non-degenerate closed two-form, i.e., a symplectic form. Thus, as expected, since it is 
a coadjoint orbit, the Hilbert Lie group L(5'') carries an invariant symplectic form whose value at the identity element 
1 (the constant loop equal to 1) is given by (flST l. 

Now we introduce the Hilbert transform on the circle 
-+ n.?«(.?lcot = h 



,Jifu(e):= — f dsu(s)cot[ = — -f diM(0 -i)cot - := lim - / dsu(e - s)cot( -] (16) 

2kJ-ti \ 2 J 2kJ-k \2/ £^0+ K Je<\s\<K ^2/ 



for any u G L^{S^), where / denotes the Cauchy principal value. We adopt here the sign conventions in I Kingl 20091 
Formulas (3.202) a nd (6 38), V ol. 1]. If m S L^iS^), then Jffu G L^iS^) and it is defined for almost every G [-7r,n:] 
(Lusin's Theorem, iKingl 120091 §6. 19, Vol. 1]). The Hilbert transform has the following remarkable properties that will 
be used later on: 

• If m(0) = i:~^_„M(n)e'"^ eL^(5'), where m(«) := ^ f^^dB u{B)e-'"^ , so = m(-«) since ii is real valued, 
then 
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Jffu{e) = -i ^ M(n) sign(«) e'"^ eL^ (5'^ 



(17) 



which follows from the identity ^/(«) sign(«) ( rtKind,l2009l Formulas (6.100) or (6.124), Vol. 1]). Here, 

sign(«) = 1 if « £ N, sign(«) = — 1 if « G — N, and s ign(O) = 0. N ote that J^m is also real valued since u{n) sign(«) = 
—u{—n) sign(— «). The formula above implies that (I King . 20091 Formula (6.126), Vol. 1]) 

• For every u G Lr{S^), we have the orthogonality property (I Kind 20091 Formula (6.127), Vol. 1]): 
Take the orthonormal Hilbert basis {fn{0) ■= e'"^ \ neZ} of L-{S^). Then ( iKingll2009l Formula (6.131), Vol. 



1]): 



Jf(Pnie) ^ -isign{n)(p„{e), for all « G Z. 



So, the eigenvalues of J^f are: — i for a ll n > 0, i for all n < 0, and if « = 0. 
• If M,v G L^iS^) then ( iKingll2009l Formula (6.99), Vol. 1]) 

(m,v) = ^ diM(i)^ d.?v(i)^ +(^M,Jfv) 

and hence ( lKingll2009l Formula (6.97), Vol. 1]) 



Il«lli2(5.) = J_^dsu{s) j + ||^«ll22(5,) 

for any u G L^(5''). This shows that || J^h||^2(51) ^ ll"ll^2(5i) the constant 1 is the best possible (I Kingl l2009l 

Formulas (6.167) and (6.168), Vol. 1]). In particular, if the average of u is zero, then is an isometry of L^i S^). 

• The Hilbert transform is skew-adjoint relative to the (5' )-inner product, i.e., J^* = —J^/f ( [Kingj2009l Formula 
(6.98) or (6.106), Vol. 1]). 

• For any u G L^{S^) we have ( llKingll2009l Formula (6.34), (6.82), or (6.156), Vol. 1]); 

Jif\(e) = -m(0) + -J- r dsu(s) = -m(0) + m(O). 
2n J-n 

• For any u G H'{S^) with s > we have M'u G H'iS^); this is an immediate consequence of ( fTTI ). If s > 1, then 



M'u' ^{,yfu)', i.e., Jf'c 



cie ~ cie 



)^on//"(5') withs> 1. 



Using these properties, if m(0) = Y.n=-ooU{n)e'"^ G H\S^), then u'{e) = Y.T,=-^u{yi)me'"^ G and hence 



(jfM')(0) = (^M)'(0) 

On the other hand, if v G ), then 



^ M(n)sign(n)e" 



and hence if m G //' (5' ^ 



^v(0)= n^v{n)e 



^ \n\u(ri)e" 



(18) 



(19) 



(20) 
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by (fTsl l. By the previous properties we have o d / dd)^ = —dP-jdO^, as expected; note that the the extra term, which 
is the zero order Fourier coefficient, does not appear in this case, because the derivative eUminates it. 

Now, if (p = e'-' e L(5'), i.e., (p(l) = 1 and / : [-tt, tt] — 5- R is a periodic function, then /(O) = /(O) = 0. Similarly, 
if M G L(M), i.e., m(1) = and we think of u as a periodic function u : [—n,n] — > K, then u{0) = m(0) = 0. This, and 
the properties of the Hilbert transform on the circle, imply: (L(M)) C L(R), is unitary on L(R) (relative to the 
//'-inner product), o = — / on L(R). Concretely, the Hilbert transform on L(R) has the form: 



m(0) 



;iGZ\{0} 



,y^u{e) = -i ^ i?(«)sign(«)e'"'' e L( 

«gz\{0} 



Thus, defines the structure of a complex Hilbert space on L(M), relative to the //* inner product, i > 1. Hence, 
translating to any tangent space of L(5' ), we obtain an invariant almost complex structure on the Hilbert Lie group 
L(5') which i s, in fact, a co mplex structure. For general criteria how to obtain complex structures on real Banach 
manifolds, see lBeltita 1 2005 1 : the argument above is a very special case of these general methods. 

Finally, L(5i) is a Kahler manifold, as proved in lAtivah and PresslevI 11198311 . This is immediately seen by noting 
that 



g{\){u,v):^(0{y^u,v)^ ^ |«|;?(«)v(«) 



(21) 



is symmetric and positive definite and so, by translations, defines a weak Riemannian metric on L(5''). Note that this 

metric is not the //' metric for any s> 1. In fact, the metric g is incomplete, whereas the //■' metric is complete . 

Concluding, (L(5M,0,g, J^) is a weak Kah l er man ifold and all structures are group invariant (see lPresslev 1 1982 1. 
Ativah and PresslevI il983n . |Presslev and Segall il986ll ). 



3.4 Weak Riemannian metrics on L(5^) 

Th e three metrics discussed in fj2]for L(5''), viewed as a coadjoint orbit of its central extension, have been computed 



by lPressleyl 11198211 . We recall here relevant formulas. 

The induced metric is defined by the natural inner product on L(R), which is the usual L^-inner product. Hence, the 
induced metric is obtained by left (equivalently, right) translation of the inner product 

b{\){u,v) :- (m,v) = ^ j\tu{t)v{t) (22) 

for any two functions m, v £ L(M). 

Define the following inner products on L(M): 

b2{l){u,v):=b{l){u,.3^v')^{u,M'v'), if u,v eH'iS^), s>\ (23) 
bi{l){u,v):=b{l){u',v')^{u',v'), if u,v eW{S^), s>l. (24) 

BiUnearity and symmetry of b\{\) and b2{l) are obvious. If m G L(5''), writing u{d) = T,n=-ooU{n)e^"^ with m(0) — 0, 
we have u'{9) = i^~^_„ni?(«)e'"^. Since {e'"^ | « G Z} is an orthonormal Hilbert basis of L?{S^), we get 

bi{l){u,u)= «V(n)P>0. 

;7 = — oo 

In addition, bi{l){u,u) = if and only if u{n) ~ for all « ^ 0, i.e., u{9) — m(0) = 0. This shows that ^^i (1) is indeed 
an inner product on L(R) which coincides with the //' inner product. Hence, if L(R) is endowed with the //^ topology 
for i > 1 , this inner product is strong if = 1 and weak if s > 1 . Left translating this inner product to any tangent space 
of L(5'') (endowed with the topology for i > 1), yields a Riemannian metric on L(5'') that is strong for 5=1 and 
weak for s > I. This Riemannian metric is the normal metric on L(5''). 

The inner product ^2 ( 1 ) is identical to g( 1 ) by ( I2TI ). ( |23l l, and ( fTSl ). Thus, translating this inner product to the tangent 
space at every point of the Hilbert Lie group L(5'' ), yields the standard Kahler metric b2 = g on L{S^ ), endowed with 
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the W topology for 5 > 1. Note that if m G L(5'), then 

which shows that the Kahler metric b2 coincides with the //'/^ metric and is, therefore, a weak metric on L(5''). 
There are relations similar to (fTSI l and (fT4l) . namely 

b{\){u,v)^bi{l){s^\,v), ^2(l)(«,v) =Z7i(1)Km,v), 

where 

/; = — 00 /7 = — 00 

if u{0) = ^~^_„M(«)e'"^. However, note that the relation involving jz/^ requires that u G H''{S^) with i > 2. 
3.5 Vector fields on L(5') anrf L(M) 



Recall that the exponential map exp : L(]R) 9 m H> e'" G L(5'') is a Lie group isomorphism ( iPresslev and Segal , [l986[ 



page 151, §8.9]). Here, we identified the Lie algebra of with M, even though, naturally, it is the imaginary axis, the 
tangent space at 1 G 5' to 5' . This means that care must be taken when carrying out standard Lie group operations with 
the exponential map, interpreted as the exponential of a purely imaginary number Since such computations affect our 
next results, we clarify these statements below. 

The tangent space at the identity 1 to is the imaginary axis. This is the natural Lie algebra of the Lie group 
and the exponential map is given by exp : iR 9 (be) 1— ^ e"^ G 5^ Of course, traditionally, one identifies ilR with M by 
dividing by i and thinks of the exponential map as exp : M 9 x g 5^ Unfortunately, this induces some problems. 
For example, since (left) translation is given by LgLie" :~ e^e", it follows that 



riL,.(iy):=-| 



de 



iye^\ (25) 

£=0 



£=0 

SO the identification of the Lie algebra with R poses no problems and we have, dividing both sides by i, 

TiL^..{y)=ye'\ (26) 
However, the definition of the exponential map for any Lie group G with Lie algebra g, yields 

^exp(f^) = r«Lexp(,^)^, forall ^ G fl. (27) 

This formula works perfectly well if the Lie algebra of 5' is iR. Indeed 

d 

— t 
dt 

which coincides with ( |27] i in view of ( |25] l. On the other hand, if the Lie algebra is thought of as R, i.e., the right hand 
side needs to be divided by i, then with the definition of exp(f.ic) = e"'' the identity above is no longer valid. What we 
should get is 

— exp(fx) = xexp(fx) = TiL^^^(,^)X = xe"^ 

by ( l26b if exp(fx) = e'", but the right hand side gives ixe"^, as we saw above. In other words, if the Lie algebra of is 
thought of as R, as is traditionally done, then we need a formula for the derivative of the Lie group exponential map in 
terms of the exponential map of purely imaginary numbers. In view of the previous discussion, this formula is 
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dt ^ ' 1 dt 

With these remarks in mind, we shall now compute the push-forward of a vector field on L(M) to L(5''). 
Proposition 2. Let X G X(L(K)) be an arbitrary vector field . Then its push-forward to L(5')) has the expression 

(exp,X)(e'") =Xiu)e"' 

for any u G L(M). 

Proof. By the definition of push forward of vector fields by a diffeomorphism, we have 



(exp.X) [e"'j = (rexpoXoexp-i) [e'"j = T,exp{X{u)) = — 



£=0 



d 

'-X{u)e' 



exp(M)exp(eX(M)) = ( — 



exp(eX(M)) ) exp(M) 



exp(M + eX(M)) 

(ID / 1 d 



£=0 



i de 



£=0 



as stated. 



3.6 The gradient vector fields in the three metrics ofh{S 



We compute now the gradients of a specific function using the three metrics. 
Theorem 1. The gradients of the smooth function H : L(5' ) — > M given by 



He' 



1 

47r 



d0/'(0) 



(i) V^H (e'-^) = fe^f for the normal metric bi; 

(ii) WH (e'-^) = — /"e'-^ with respect to the induced metric b for f G W{S^ ) with s > 2; 

(iii) V^i/ (^'■^) ~ with respect to the weak Kahler metric b2- 

Proof, (i) Since T\L^ifU ~ ue'f for any u G L(M) and e'^ G L(5''), invariance of bi yields 



bi{l) (e-''^V^H(e'A ,m) ^b 



d 
dt 



dH e 



t=0 



i- Tde (/'(0)+f«'(0))' 

f=0 — 



which shows that V^// (e'-^^) = /e'-/^. 

(ii) Proceeding as above, using the same notations, and assuming that/ G//" (5') withi>2, we have 



^(1) (e-'-^V// (e'f^ ,i^=b {e'f^ (VH 



e'f] .ue'f 



dH ( e'n ( ue'f 



= (-/",«)®Ki)(-r,«) 



which shows that V// (e'/) -/"e*'''. 



12 Anthony M. Bloch, Philip J. Morrison, and Tudor S. Ratiu 

(iii) This computation uses the isometry property of relative to the iJ inner product. We have, 
biiX) [e-'fy^H [e'f^ ,uj=b2 (e'^) (v^// (e^/) = dH (e'^] ( ue'f^ 



which shows that V^H (e'f) = {Jff')e'f. 
Since 



(0 ( e'f) ( Ji-V^H (e'f) ,ue'f) ^ b2 (e^A ( V^// (e'f) ^f) = dH (e^A (ue'f 



it follows that the Hamiltonian vector field on (L(5''),a)) for the function H is Xh = J^V^H. Since commutes 
with the tangent lift to group translations, TheoremlTJiii) implies that 



This proves the first part of the following statement. 

Corollary 1. The Hamiltonian vector field ofH relative to the translation invariant symplecticform CO on L(5' ) whose 
value at the identity element is given by dlSI l has the expression Xh (e''^) = —fe^'^- Its flow is the rotation 



Proof. Since L(]R) 3 u \ — > e'" G L(5' ) is the exponential map and we think of R as the Lie algebra of 5' (and not the 
imaginary axis), we write de^'" /dt = we"" without the factor of i in front (see (|28])). The verification that Ft is indeed 
the flow of Xh is straightforward: 



= XH(FJe^f)){e) 



as required. 



We recover thus IPresslevll 19821 Proposition 3.1] (up to a sign which is due to different conventions calibrating (O, 
Jif, and b2). 

Applying Proposition [2] to TheoremfT] we get the following result: 
Corollary 2. The three gradient vector fields for the smooth function Hi : L(R) R given by 

Hi{u) = ^pAe{u'f 

are 

(i) V'//i (u) — ufor the weak inner product bi{\) defining the normal metric; 

(ii) V//i (m) = ~u" for the weak inner product b{l) defining the induced metric, where for u G //''(M) with s > 2; 

(iii) V^Hi (m) — Jffu' for the weak inner product ^2(1) defining the Kdhler metric. 

Since the exponential map is a Lie group isomorphism and the three metrics coincide with the respective inner 
products at the identity, their left invariance guarantees that the three inner products on L(M) correspond to the three 
invariant metrics on L(5'). 

Applying Proposition|2]to Corollary[Tl we conclude: 

Corollary 3. The Hamiltonian vector field of Hi relative to the symplectic form CO given by (115b has the expression 
Xh{u) = -m'. Its fiow is {F,(u)) (0) = u{e-t). 
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The verification of the statement about the flow is immediate: 
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^ (Friu)) (9) = ju{e -t) = -u\e -t) = {Xh [Friu))) (0). 

If one is willing to put more stringent hypotheses on the functional, it is possible to obtain a general result. 

Theorem 2. Let H : L(5') -> R fee a smooth function {with L(5') endowed, as usual, with the //' topology for s>\) 
and assume that the functional derivative 5H / 5u G L(5') exists. Then the gradient vector fields are 

(i) VH{u) = ^ with respect the weak inner product b(\) defining the induced metric; 



(ii) (v1//(m)) (0) = -HqAcp (^/(f dv/ff with respect to the (weak) inner product b\{\) defining the normal 
metric, provided both jQd(p^{(p) and J^dq) (^J^dxj/ are periodic; 

(iii) {V^H{u)) (0) = -jr/(fd(p|^((p) with respect to the weak inner product ^2(1) defining the Kahler metric, 
provided j^dip ^ ((p) is periodic. 

Proof, (i) For the inner product b{\) on L(5'') defining the induced metric, if m, v G L(M), we have by periodicity of 



„V„(.,„.).D«,„),..(|^,.)raK.)(^,.) 



This shows that V//(m) = 

(ii) For the inner product b\{\) on L(5') defining the normal metric, if «,y £ L(R), we have by periodicity of 



i„(l)(V«(„),,.) = DH(.).v= ('I^.A = 



de 



which shows that {V^H{u)){B) = -J^d(p (^/g^'dv/ff (r)^ 

(iii) For the inner product ^2(1) on L(5'') defining the Kahler metric, if u,v G L(R), we have by periodicity of 
[(fdip ^{(p) and the isometry property of , 



fe2(l)(V"//(M),v)=D//(M). 



\ 8u 



d0— (0M0) 

2n J-n OU 



OU 



1 

. 2n J-n 



6 8H 
d<p (<p) ) v'(0) 

OM 



/'d(p^((p),^v'\ 1*^2(1) f-^^ /'d(p^((p),V 

OM / V Jo OM 



which shows that {V^H{u)) (9) = -.Jf £d(p ^{(p). U 

Corollary 4. Under the same hypothesis as in Theorem |2jiii), the Hamiltonian vector field of the smooth function 
H ; L(5') — > R relative to the symplectic form CO on L(R) given by il5i has the expression Xh{u) ~ /o^d(p ^(<p) 



Proof We have Xniu) = Jify^H{u) /(f d^ |f (^). 
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Of course, using Proposition |2] there are immediate counterparts of Theorem |2] and Corollary |4] on the loop group 
L(5''), which we shall not spell out explicitly. 

The hypotheses guaranteeing the existence of the functional derivative of H relative to the weakly non-degenerate 
pairing are quite severe. For example, the theorem can be applied to the functional Hi in Corollary[2] but one needs 

additional smoothness. Indeed, the first thing to check is if this functional has a functional derivative. In fact, it does 

not, unless we assume that u G H\S'^) for i > 2, in which case we have 



D//i(„) .v = -L r dsii'{s)v'{s) = ^u'{s)v{s) 
2n J -71 2n 



1 



i.e., 5H / 5u = — m". With this additional hypothesis, the gradient flow with respect to the weak inner product b{l) 
defining the induced metric is given by «, ~ —u". 

Therefore, to continue computing the other two gradients of H\, we need to assume that u e H'{S^) for s>2. 
Provided this holds, to find the gradient relative to the (weak) inner product ( 1 ) defining the normal metric, we have 
to check that both 

8h 

/ d(p^{(p) = ~ d(pu"{(p) = -u'{e) + u'{0) 
Jo on Jo 

f-e / r(p 5// \ f-e 

d(p n d\if — {^) \ = -J^ d(p{u'{(p)~u'{0))^-ui9) + u'{0)G 

are periodic. While the first one is periodic, the second one is not unless we assume that m'(0) = 0. With this additional 
hypothesis, the gradient is given by m, = u. However, we know from Corollary |2] that neither i > 2, nor m'(0) = is 
needed. In addition, this can also be seen directly, as follows. For any u,v e L(M), we have 

/7i(1)(V1//(m),v) =D//(m)-v= r dsu{s)v'{s) = (u' y) ^ bi{u,v) 

2u J 

which shows that V'//(m) = u. 

The same situation occurs in the computation of the third gradient. In the hypotheses of the theorem, we have 

{V^H{u)){e) = -J^ / d(p—{(p)^.J'/f(u'-u'{0))^Jfu' 



8 II 

because the Hilbert transform of a constant is zero. Thus, the gradient flow is given in this case by 

Uf = '^u' [ y 



As before, the same result can be obtained easier and without any additional hypotheses in the following way; 

b2{l){y^H{u),v) = DH{u) - v ^ {u ,v') = {jifu ,M'v')^ b2{l){J^fu' ,v). 



3.7 Symplectic structure on periodic functions 

The form of the periodic Korteweg-de Vries (KdV) equation we shall use is 

ut-6uue + ueee=0, (29) 

where u{t, 0) is a real valued function of f G M and 9 G [— tt, n], periodic in 9, and ug :~ du/d9. The KdV equation 
is, of course, a famous integrable infinite dimensional Hamiltonian system. It is Hamiltonian on the Poisson manifold 
of all periodic functions relative to the iGardner. Ii,1971il bracket 
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= — / d0— — — , (30) 
2k J-ji ou dO ou 

where 

F{ii) = / d0/(M,Me,Mee,...) 
and similarly for G; the functional derivative 5F / 5u is the usual one relative to the L?'{S^) inner product, i.e., 



5F 






fdf_\ 


d^ 


( df \ 


5u 


du 


d9 


ydue J 


+ d¥ 





The Hamiltonian vector field of H{u) = [^^,.(19 h{u,ug,ugg, . . .) has the expression 

d f5H\ 



For the KdV equation one takes 



da \ ou J 

The Casimir functions of the Gardner bracket are all smooth functionals C for which 5C/5u = c is a constant 
function, i.e., 

1 



C(u)=^(c,u) = — / decu{e)=cu(0). 

2k J -71 

Thus C-\Q) is a candidate weak symplectic leaf in the phase space of all periodic functions. The situation in infinite 
dimensions is not as clear as in finite dimensions, where this would be a conclusion, because there is no general 
stratification theorem and one cannot expect, in general, more than a weak symplectic form. However, in our case, this 
actually holds, as shown in Zaharov and FaddeevI 1 1971 1. Indeed, 

a{uuu2) ' = -^1 (^Jq'^^ (mi(<p)m2(0)-m2(<p)mi(0))^ ""2^/ i^Io '^^"i('P)) "2(0) 

d(j!)Mi((j!)),M2^ (32) 

defines a weak symplectic form on L(IR) whose formal Poisson bracket is ( l30b . This immediately shows that there is 
a tight relationship with the symplectic form (O of the complex Hilbert space L(M), the Lie algebra of the based loop 
groups, given by ( fTSl ), namely 

n 

for all M, V G L(M) of class H\ s>2. Defining 



-j «:=y^d<p«(<p), 
the KdV symplectic form a has the suggestive expression (see (l28l l) 



which is well defined onH ^ (5' , M) . 



On the other hand, the Poisson bracket given by the Kahler symplectic form (flSl l on L(R) is 



1 ,„5F f d 5G 



^'^''^-2-KL''j-u\Te J-u' '''' 
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which 
gli.e.. 



which is similarly well defined onH ^, and the Hamiltonian vector field defined by this bracket is given by Corollary 



Now, the gradient vector field for the corresponding Kahler metric, as computed in Theorem|2|iii), is written as 

, f d\-^ 5H 



4 Metriplectic Systems 

In this section we define metriplectic systems and show how to construct general classes of such systems in terms of 
triple brackets for both finite- and infinite-dimensional theories. We use some of the machinery developed above to 
address specific examples. 



4.1 Definition and consequences 



A metriplectic system consists of a smooth manifold P, two smooth vector bundle maps tt, k : T*P TP covering the 
identity, and two functions H,S <E C°°{P), the Hamiltonian or total energy and the entropy of the system, such that 

(i) {F, G} := {dF, n{dG)) is a Poisson bracket; in particular n* ~ —n; 

(ii) {F,G) := {dF,K{dG)) is a positive semidefinite symmetric bracket, i.e., (,) is R-bilinear and symmetric, so 
K* = K, and {F,F) > for every F e C~(P); 

(iii) {S,F} = and {H,F) = for all F e C°°(P) <^ 7t{dS) = K{dH) = 0. 



The metriplectic dynamics of the system is given in terms of the two brackets by 

jF = {F,H + S} + {F,H + S)^{F,H} + {F,S), for all FeC~(P), 
or, equivalently, as an ordinary differential equation, by 



(36) 



-c(0 = n{c{t))dH{c{t)) + K{c{t))dS{c{t)). 



(37) 



The Hamiltonian vector field Xh n{dH) G X(P) represents the conservative or Hamiltonian part, whereas 
K"(d5') e X.{P) the dissipative part of the full metriplectic dynamics ( l36b or dSTl i. 



As far as we know , first atte mpts to introduce such a struct u re were given in adjacent papers by iKaufmanI 11198411 
and lMorrisonI ll 1 984all . (See also Kaufman and Morrison 1 1982 1.') Kaufmanl l 1984 1 imposed, instead of (iii), the weaker 
condition {H,S} = {H^S) = 0, which is enough, as will become apparent below, to deduce the First and Second Laws 
of Thermodynamics. In the plasma examples presented, he used (iii) for a large class of fun ctions. All t hree ax ioms, 
including the degeneracy condition of (iii), were stated expU citly in iMorrisonI il984all and iMorrisonI lll984bll . The 
former treated the same kinetic example as iKaufmanI Ml 98411 along with additional formalism, while the latter pre- 
sented the metriplectic formalism for th e compressible Navier-Stokes equations with entropy production. All three 
axioms were restated in Morrison 1 1986ll . where the terminology metriplectic was introduced and a detailed physical 
motivation for the introduction of (iii) is presented along with other examples such as a dissipative free rigid body 
equatio n and the Vlasov-Poisson equ ation with a collision term that generalizes the Landau and Balescu-Lenard equa- 
tions. In lGrmela and Ottingen 1199711 . under the name GENERIC (General Equations for Non-Equilibrium Reversible 
Irreversible Coupling), the same geometric structure was used to analyze many other equations; due to this paper and 
subsequent work of these aut hors, the metrip lectic formalism has been popularized. For a very interesting modern ap- 
pUcation of this structure see iMielkd 1201 ill and for further discussion about avenues for generalization see iMorrison 
1 2009ll . 
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The definition of metriplectic systems has three immediate important consequences. Let c(f ) be an integral curve of 
the system (|37] |. 



(1) Energy conservation: 

(2) Entropy production: 



dt 



H{c{t)) = {H,H}{c{t)) + {H,S){c{t)) = 0. 



dt 



Sic{t))^{S,H}{c{t)) + {S,S){c{t))>0. 



(38) 



(39) 



(3) Maximum entropy principle yields equilibria: Suppose that there are n functions Ci , . . . ,C„ G C°°{P) such that 
{F,C,} = {F,Ci) = for all F g C°°{P), i.e., these functions are simultaneously conserved by the conservative and 
dissipative part of the metriplectic dynamics. Let pQ G P be a maximum of the entropy S subject to the constraints 



H-^{h) n Cf ' (ci ) n . . . C,7 ' (c„), for given regular values /i , ci , . . . , c„ G M of //, Ci , . . . , C„, respectively. By the La- 
grange Multiplier Theorem, there exist a, jSi , . . . , j3„ G M such that 

dSipo) = adH{po) + jSidCi (po) + ■ ■ ■ + dC„{po). 

But then, assuming that a ^Q, for every F G C°°(P), we have 

{F,H}{po) + {F,S){po) = {dFipo),7z{po) (d//(po))) + {dF{po).K{po) {dS{po))) 



= ( dF{po), -7t{po) (dSipo) - j3idCi (po) dC„{po)) 

- (dF(po), K{po) (ad//(po) + jSidCi (po) + • ■ ■ + dC„(po))> 



l^{m{po)-^{F,C,}{po)- 
+ aiF,H){po) + lii{F,Ci)ipo) 



■ + l3„iF,C„){pQ) = 



which means that po is an equilibrium of the metriplectic dyna mics (l36ll or ([37ll. Th is is akin to the free energy 
extremization of thermodynamics, as noted by IMorrisonI il984bll and lMorrisonl il986ll where it was suggested that 
one can build in degeneracies associated with Hamiltonian "dynamical constraints." (See also|Mielke |2011].) 

Suppose that K ^C°° (P) is a conserved quantity for the Hamiltonian part of the metriplectic dynamics, i.e., {K,H} — 
0. Then, if c{t) is an integral curve of the metriplectic dynamics, we have 



dt 



K{cit)) = dK{c{t)) {c(t)) = {dF{c{t)), n{c{t)) (d//(c(f )))) + {dF{c{t))Mc{t)) (d5(c(f )))) 
= {K,H}icit)) + iK,S)icit))^iK,S)icit)). 



As pointed out in IMorrisonI ill 986ll . this immediately implies that a function that is simultaneously conserved for the 
full metriplectic dynamics and its Hamiltonian part, is necessarily conserved for the dissipative part. Physically, it 
is advantageous for general metriplectic syst ems to conserve dynamical constrain t s, i.e., con s erved q uantit ates of its 
Hamiltonian part and the examples given in KaufmanI 1 1984 1. Morrison 1 1984a l. Morrison 1 1984bll . and MorrisonI 
lfl986ll satisfy this condition. 



4.2 Metriplectic systems based on Lie algebra triple brackets 

Associated with any quadratic Lie algebra (i.e., a Lie algebra admitting a bilinear symmetric invariant form) is a natural 
completely antisymmetric triple bracket. This is used to construct Lie algebra based metriplectic systems. The algebra 
so (3) is worked out explicitly and examples are given. 
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4.2.1 General theory 

A quadratic Lie algebra is, by definition, a Lie algebra admitting a bilinear symmetric non-degenerate invariant form 
ff : X g — > R (the letter K is meant to remind one of the Killing form in a semisimple Lie algebra). Recall that invariance 
means that K'([^ , 77], Q = K{£,,[rj,Q) for all ^,J7,C G or, equivalently, that the adjoint operators adtj for all 77 G g 
are antisymmetric relative to K. Non-degeneracy (strong) means that the map g 9 ^ 1-^ k(£, , ■) € q* i s an iso morphism. 
Finite dimensional quadratic Lie algebras have been completely classified in [Medina and Revovl 11 198511 . For finite 
dimensional Lie algebras, non-degeneracy is equivalent to the following statement: K{E,,ri) =0 for all 7] £ g if and 
only if ^ = 0. In infinite dimensions this condition is called weak non-degeneracy and it is implied by non-degeneracy 
but the converse is, in general, false. 

For example, let g be an arbitrary finite dimensional Lie algebra. Recall that the Killing form is defined by K'(^ , 77 ) : = 
Trace(ad^ o ad^j ). If {e,}, i—l,... dimg, is an arbitrary basis of g and c'^^j are the structure constants of g, i.e., [e,, e^] = 
c'^ijSp, then 

and hence the components of k in the basis {e,}, / = 1 , . . . dimg, are given by 

Kij = K{ei,ej)^cP.^c''jp. 

The Killing form is bilinear symmetric and invariant; it is non-degenerate if and only if g is semisimple. Moreover, — K 
is a positive definite inner product if and only if the Lie algebra g is compact (i.e., it is the Lie algebra of a compact Lie 
group). 

In general, let K" be a bilinear symmetric non-degenerate invariant form and define the completely antisymmetric 
covariant 3 -tensor 

In the coordinates given by the basis {e,}, / = 1, . . .dimg, the components of c are 



This construction immediately leads to the triple bracket introduced bv lBialvnicki-Birula and MorrisonI 11199 111 (see 
also lMorrisonI il998ll ). {•,•,•}: C°°(g) x C°°(g) x C"(g) C"(g) defined by 

{f,g,hm--=c{VM),Vg{^),Vh{^)):^K{Vm,[Wg{^)yhm, (40) 

where the gradient is taken relative to the non-degenerate bilinear form K, i.e., for any ^ G g we have 

or, in coordinates 

where [k'^] ~ [k^/]^', i.e., K'^Kji^ = Sj^. This triple bracket is trilinear over M, completely antisymmetric, and satisfies 
the Leibniz rule in any of its variables. In coordinates it is given by 

{f,g,h} = C,,,V'fV^gV'h = K,,,c'\K''^^K^''^K'rdh ^ jc^krdfdgdh 

u>6, J ijK J s jk ^^,1 jk d^P d^i d^'' 

nc,r dg dh 



where c^'''' are the components of the contravariant completely antisymmetric 3-tensor c associated to c by raising its 
indices with the non-degenerate symmetric bilinear form K, i.e., for any ^ , Tj, 1^ S g, we have 



c(K-(.^,-),K-(r7,-),K-(7r)) ■■=c{^,ri,Q. 
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This construction extends the bracket due to iNambul 11197311 to a Lie algebra setting. Nambu considered ordinary 
vectors in and defined 

[f.gMmmbuin) = v/(ir) • (v^(n) x vh{n)) , (4i) 

where ' •' and ' x ' are the ordinary dot and cross products. Thus, the Nambu bracket is a special case of the triple bracket 
( |40] i in the case of g = so(3), whose the structure constant s are the completely antisymmetric Levi-Civita symbo l Sijic- 
Such 'modified r igid body brackets' were also described in lBloch and Marsdenljl 99011 . iHokn and Marsdenllll99lll . and 



Marsden and Rat iu [19993. 

If g is an arbitrary quadratic Lie algebra with bihnear symmetric non-degenerate invariant form K, the quadratic 
function 

C2{^):=\Ki^,^) (42) 
is a Casimir function for the Lie-Poisson bracket on g, identified with g* via K, i.e., 

{f,g}±i^)^±K{^,[vm,vgm, (43) 

as an easy verification shows since VC^i^) = ^. In view of ( l43T l, the following identity is obvious 

{/,§}+ ^{C2,f,8} 

(this was first pointed out in Bialynicki-Birula and MorrisonI lll99lll ). For example, if g ~ so (3), the (-)Lie-Poisson 
bracket 

{f,gY_"^'\n) = -{C2,/,g}Nambu(iT) = -H • (v/(n) X vgin)) (44) 

is the rigid body bracket, i.e., if h{n) — jFI £2, where 77, = /,i2,, /, > 0, / = 1,2,3, and /, are the principal moments 

of inertia of the body, then Hamilton's equations ^F{n) ~ {f,hY_°^^\n) are equivalent to Euler's equations fl = 
nxQ. 

Note that given any two functions, f,gG C°°(g), because the triple bracket satisfies the Leibniz identity in every 
factor, the map C°°(g) B h i—> {h,f,g} G C°°(g) is a derivation and hence defines a vector field on g, denoted by 
Xf,g : g ^ g, i.e., 

{dhi^),Xf,g{^)) = K{Vh{^),Xf^^i^)) ^ {h,f,gm for all /i GC~(g). (45) 

Note that Xfj = 0. Thus, for triple brackets, two functions define a vector field, analogous to the Hamiltonian vector 
field defined by a single function associated to a standard Poisson bracket. 
From ( l40l i we have the following result. 

Proposition 3. The vector field Xfg on g corresponding to the pair offiinctions f,g is given by 

Xf.,{^)^[Vm,Vg{^)]. (46) 

Triple brackets of the form ( |40l i can be used to construct metriplectic systems on a quadratic Lie algebra g in the 
following manner. Let k be the bilinear symmetric non-degenerate form on g defining the quadratic structure and fix 
some h G C°°(g). Define the symmetric bracket 

[f,g)l{^):^-K{X,j{^),X,,,{^)). (47) 

Assume that — JC is a positive definite inner product. Then (/,/) > 0. Thus we have the manifold g endowed 
with the Lie-Poisson bracket ( |43] |. the symmetric bracket (l47T i. the Hamiltonian h, and for the entropy S we take 
any Casimir function of the Lie-Poisson bracket. Then the conditions (i)-(iii) of i i4.1l are all satisfied, because 
{h,g)^ = —k{Xi, i,,Xi, g) = -~K{Q,Xi,g) ~ for any g G C°°(g). The equations of motion ( |36] | are in this case given 
by 

j/i^) = K (v/(^), 1^) = {fM±{^)+ if,sm = ±K{^,[vmyhm -K{x,ji^),xi,sm 

= T'c(V/(^),[^,VM^)])->c([V/,(^),V/(^)],[VM^),V5(^)]) 
forany/GC"(g). 
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This gives the equations of motion 

^ =±[^,VK^)] + [VK^),[VK^),V5(^)]]. (48) 

Note that the flow corresponding to 5 is a generalized double bracket flow. Observe also that this flow reduces to a 
double bracket flow and is tangent to an orbit of the group if V/i(^ ) = ^ . Indeed if /j = j , ^ ) the symmetric bracket 
WT\ reduces to the symmetric bracket induced from the normal metric. 



4.2.2 Special case of so (3) 

If the quadratic Lie algebra is so(3), we identify it with M? with the cross product as Lie bracket via the Lie algebra 



isomorphism ' 



so (3) given by uv := u x v for all u, v G R^. Since Ad^u = Au, for any A E 5(9(3) and u G 



we conclude that the usual inner product on M? is an invariant inner product. In terms of elements of so (3 ) we have 
u- V = — J Trace (uv). We shall show below that the metriplectic structure on Sr" is precisely the one given in lMorrison 

lim. 



Recall that the Nambu bracket is given for so(3) by (HD) and hence the symmetric bracket (l47b has the form 



K{{n,ii,f},{n,h,g})^e'' 



dh df dh dg 

'J 



Jnm St 

'■ dn" dn" dw 

\VhfVg-Vf-{Vf.Vh){Vg-Vh) 



(49) 



where in the third equality we have used the identity e"""£.^' = _ This coincides with I Morrison , 19861 

equation (31)]. 

With the choice S{n) — |n |P/2 and the usu al rigid body Hamiltonian, the equations of motion ( |48] | are those for 
the relaxing rigid body given in lMorrisonHl986ll . 



Comments. 

• In three dimensions any Poisson bracket can be written as 



dW dW 



df dg 

dn' dw 



(50) 



where = 1,2,3, and V G . The last equality follows from the identification of 3 x 3 antisymmetric matrices 
with vectors (the hat map discussed above). Using the well know fact (which is easy to show directly) that brackets 
of the form of (ISOl l satisfy the Jacobi identity if 



we conclude that 



y-Vxv = 0, 



{F,G}f — {/,/^,G}Nambu 



(51) 

(52) 



satisfies the Jacobi identity for any smooth function/; i.e., unlike the general case where the theorem of lBialvnicki-Birula and Morrisoi 
1 199111 requires / to be the quadratic Casimir, one obtains a good Poisson bracket for any /. Thus, for the special 
case of three dimensions, one can interchange the roles of Hamiltonian and entropy in the metriplectic formalism. 



Thinking in terms of so (3)*, the setting arising from reduction (see e.g. Marsden and Ratiu 1 1999ll X this construction 
leads to a natural geometric interpretation of a metriplectic system on the manifold f = . With the Poisson bracket 
on M? of (|52l), the bundle map n : T*R^ TR^ has the expression 



;r/(x,n) = (x,V/(n)x(.)T 
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since d//(iT)^ = VH{n) (d//(i7) is a row vector and VH{n) is its transpose, a colum n vector). Now the triple 

bracke t associated to the equation (|48] | can be used to generate a symmetric bracket given in lBloch. Krishnaprasad. Marsden. and Ratiu 
I I994I1 as follows: 



{F,G)BKMR{n) = {F,G)^ = K{{n,c,F},{n,c,G}) 
= (ii X VF(n)) • (n X VG(n)) . 

where now C = | \n\ p/2. Hence the bundle map K : T*^? TM? has the expression 

K{x,n) = -nx fiix 



(53) 



Thus, with the freedom to choose any quantity 5 = / as an entropy, with the assurance that ( ISTl ) will be satisfied 
because V x V = V x V/ = 0, we can take H^C and have {F,S} / = and {F,H) = for all F e C°°(M3). The 
equations of motion for this metriplectic system are 



n = -nx vf{n) -nx{nx vf{n)). 



(54) 



The sym metric bracket is the inner product of the tw o Ham iltonian vector fields on each concentric sphere. As dis- 
cussed in lBloch. Krishnaprasad. Marsden. and Ratiul il994ll . this symmetric bracket can be defined on any compact 
Lie algebra by taking the normal metric on e ach coadjoin t orbit. 
• The following set of equations were given in lFishI 1200511 : 



n = ys{n) X VH{n) - VH{n) x {VH{n) x vs{n)). 

Yet, this metriplectic system is identical to that obtained from ( |48] ). using ( |49] l. viz. 

ti = {n,5,//} + ic({n,//,n},{iT,//,5}) , 

Replacing i/ by g in ( l49l ) gives 

(F,G),((n)) = K{{{n),g,F},{{n),g,G}) = (v^(ir) x VF(ir)) • (v^(ir) x VG(ir)). 

Thus, the bundle map k: r*R^ TM? has the expression 

K-g(x,iT) = -V^(iT)x (vnx(.)T). 
Examples: Two special cases of the equation (|55] l are of interest. 



(55) 



(56) 



(57) 



(i) If we take H = ^ ||n|p and S = c ■ U, c & constant vector, we obtain 

ti = cxn-nx{nxc). 

(ii) If we take S = 7 ||n|p and H — c ■ FI, c a constant, we obtain 

n^Fixc — cx{cxn). 



(58) 



(59) 



The equations of motion (1581 ) is an instance of double bracket damping, where the damping is due to the normal 
metric, whereas ( |59] l gives linear damping of the sort arising in quantum systems. 
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4.3.1 The Toda lattice equation revisited 



We note that the Toda lattice equation fits into the metriplectic picture in a degenerate but interesting fashion since it 
has a dual Hamiltonian and gradient character which may be seen by writing it in the double bracket form (O. . 

It may be viewed either as the Hamiltonian part or the dissipative part of a metriplectic system with Hamiltonian 
H = J TtL^ or entropy function S = TiLN respectively with the Toda lattice equations in the corresponding form ^ 
or as discussed in Section|2l This observatio n may be extended to the Toda latt ice flow on the normal form of any 
complex semisimple Lie algebra as can be see in lBloch. Brocket!, and Ratiul lll992ll . 



4.3.2 Full Toda with dissipation 

It is possible to construct an interesting metriplectic system which incorporates the full Toda dynamics. 

We consider the again the flo w on the vector spac e of symmetric matrices t-^ = sX)m{n) but now consider the flow 
on a generic orbit as discussed in lDeift et al. 1 1199211 where it was shown that the flow is integrable. The Hamiltonian 



is again j TrL^ and the flow on full symmetric matrices is given by 

L=[7t,L,L] 



(60) 



with TZg being the projection onto the skew symmetric matrices in the lower triangular skew decomposition of a matrix. 
In this setting there are nontrivial Casimir functions of the bracket (|9]l. These are given as follows. For L an n x n 
symmetric matrix set for <k < [In] 



det(L- A)^. : 



n-2k 
r-O 



n—2k—r 



(61) 



where the subscript k denotes the matrix obtained by deleting the first k rows and the last colum ns. Then I\k{L) 
Eii^{L) / Eqi^{L) are Casimir functions of the generic orbit in sr)m(n) as shown in lDeift et al. I lll992ll . 
Thus we obtain the metriplectic systems 



L^[n,L,L] + [L,[L,Vhu\] 



(62) 



where the metric is the norm al metric on orbits of su(«) restricted to the symmetric matrices (identified with / times 
the symmetric matrices) as in lBloch. Brockett. and Ratiul 1199211 . Here H = ^ TrL^ and S — 



4.4 Metriplectic systems for pdes: metriplectic brackets and examples 

First we construct a class of metriplectic brackets based on triple brackets for infinite systems, then we consider in 
detail an example based on Gardner's bracket on S^. Lastly, we mention various generalizations. 



4.4.1 Symmetric brackets for pdes based on triple brackets 

Similar to il4.2l we can construct metriplectic flows for infinite-dimensional systems from completely antisymmetric 
triple brackets of the form 

{E^FM} = f deJ d02 / d03^;A(0i,02,03)(^%)(0i) i^'F^m) {^'GM) (63) 

75' 75' 75' 

where E, F, and G are smooth functions on 5', is a smooth function on 5' x 5' x 5' which is completely anti- 
symmetric in its arguments, so as to assure complete antisymmetry of {E,F,G}. In addition, we denote £„ : — 5E/5u, 
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etc. Let i~ 1 , 2, 3, be pseudo-differential operators. Evidently, the triple bracket of ( 1631 1 is trilinear and completely 
antisymmetric in £ , G. 

From ( |63] | and a Hamiltonian //, we construct a symmetric bracket as follows: 



{f,g)h^ J^de' J^de"{uie'),H,F}^{e',e"){uie"),H,G}, (64) 

where U{0) in ( |64] i denotes the functional 

u{e):u^ j de'u{e')d{e-e'). (65) 



We shall use this notation in subsequent expressions below. The 'metric' is assumed to be symmetric and positive 
semidefinite, i.e., the smooth function ^ : 5' x 5' M satisfies ^(0', 9") = ^{9", 9') and 

d9' [ d9"'^{9',9")f{9')f{9")>0 (66) 



for all functions / G C°°{S^). Therefore, by construction, it is clear that (|64] | satisfies the following: 

(i) {F,G)h = iG,F)H for a\\ F,G, 

(ii) (F,i/)//=OforallF, and 

(iii) If,F)h >0 for all F. 

As a special case suppose = £P for all 1,2,3; then ( l63l l becomes 



{£,F,G}= / d0i / d02 / d03^(0i,02,03)^(0i)£«^(02)^'„^(03)G„. (67) 
is' is' is' 

As a further specialization, suppose 'r^{9i ,92, 9^) is given by 

, 92, 03) = A(0i , 02) +A(02, 03) +A(03, 01 ) (68) 

where A is any antisymmetric function, i.e., 

A(0i,02) = -A(02,0i). (69) 

The form ( |68] |. assuming ( |69] l, assures complete antisymmetry of 

Finally, a particularly interesting, self-contained, case would be to suppose the A's come from some Poisson bracket, 
according to 

A(0i,02) = {f/(0l),C/(02)}. (70) 

It would be quite natural to choose the entropy, S, to be a Casimir function of this bracket and to choose this bracket 
as the Hamiltonian part of the metriplectic system with symmetric bracket given by (|64] |. We give an example of this 
construction in Sec. 14.4.21 

It is evident that one can construct a wide variety of symmetric brackets based on triple brackets. For example, 
one can choose the pseudo-differential operators from the list {J^j,d /d9 ,{d /d9)^^ ,J^^}, where J^^ is the identity 
operator, and the Hamiltonian, H, and entropy (Casimir) C could be one of the following functionals: 



Ho = 


( d9u 


(71) 








H2 = 


ld9u^/2 


(72) 


Hi = 


ld9.V2 


(73) 








HkcIV ~ 


J^d9 + 


(74) 
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In the Sec. 14.4.21 we will construct a metriplectic system based on the Gardner bracket ( l30l l of Sec. 13.71 To avoid 
complications, we choose a simple example, yet one that displays general features of a large class of 1 + 1 energy 
conserving dissipative system. 



4.4.2 Metriplectic systems based on the Gardner bracket 

For simplicity we choose J^, = J^^i for all /, and as mentioned above, we suppose A(0i,02) is generated from the 
Gardner bracket (l30l l. i.e., 

A(0i,02):={f/(0i),C/(02)}= / a'e 5(0 - 00-^5(0- 02) = 5'(0i- 02), (75) 

75' "0 

where prime denotes differentiation with respect to argument and 5'(0i — 02) is defined by 

\ d0i I d025'(0i-02)/(0i)g(02) = - /" d0i I AsS\s)f(Qx)^(Qi-s)= [ d0i/(0i)/(0i) 

= - / d0i / d025'(02 - 0i)/(0i)g(02) 

for any f,g e C°°{S^), which shows that 5'(02 - 0i) = -5'(0i - 02). With this choice for A we obtain 

■^(01 , 02, 03) = 5'(0i - 02) + 5'(02 - 03) + 5'(03 - 0i) , 

and Eq. ( |67] | becomes 

{E,F,G}^ [ d0i f d02 [ d03 [5'(0i-02) + 5'(02 - 03) + 5'(03-0l)] £„(0l)^'„(02)G„(03) 

Js' Js' Js' 

d0G„(0)^ d0F„(0)£;,(0) + d0£„(0)^ d0 G„(0)F,;(0) 

d0F„(0) ) /" d0£„(0)Gl(0). (76) 
s' / is' 

We shall construct a metriplectic system of the form 

F ^{H,F,G} + J^do' j^de"{u{e'),s,F}'^{e',e"){u{e"),s,G}, 

using the Gardner bracket ( TtTI i. 

Observe if we now set F = Hq, the Casimir for the Gardner bracket dTlT l. then, since SHq/Su = 1, we obtain 

{F,i/o,G}= / d0F„G;, (77) 

is' 

which is precisely the Gardner bracket. To see this, let us compute, for example, the integral in the third term of (|76] |. 
Changing variables i = 03 — 0i we get 

d0i / d02 / d03 5'(03-0l)£„(0l)G„(03) = - I ds [ d035'(i)£„(03-^)G„(03)= / d03 (03)G„03). 
Js' JS' JS' JS' JSl 

A similar computation shows that the first and second terms vanish. 

In order to construct the symmetric bracket in (|64l i. we need the following, computed using ( |76l ): 

{f/(0),i/,G} = - ^_d0G„(0)) Hl{e) + l^deG,,{e)H;,{e)+Q^deH,,{e)j g'„(0). (78) 

Now with the counterpart of dTSl l for the functional F with U{d'), a choice for //, and a choice for we can 
construct {F,G)h- We make the following choices: 
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//2(m)= / d0— , S{u):=Hq{u)= / d0M (79) 

^(0',0") = 5(0'-0"). (80) 
Now choose H2 from ( |79] l and insert it into ( |78] l which gives 

{c/(0),//2,G} = - ^^d0G„(0)^ u'{e) + j^AeGu{d)u\e)+SG'M (8i) 

and to construct the symmetric bracket ( |64l ), we need 

{f/(0"),//2,5} = -m'(0"). (82) 



Thus, the equations of motion are 



where 



jF = {F,Ho,H2} + iF,S)H2 



{F,S)^^_^ l^dO' l^dB" {U{B'),H2,F}'^{e',B"){U{e"),H2,S} . (83) 
This yields 

Ut-ug=Sugg + Q with Q := [ de'\ug,\'^ . (84) 



Equation ( |84] | has several interesting features. For fixed given constant S and it is a linear equation composed 
of the heat equation with a source and with the inclusion of a linear advection term. One can proceed to solve this 
equation by the usual method of constructing a temporal Green's function out of the heat kernel and expanding in a 
Fourier series. After such a solution is constructed, one must enforce the fact that the global quantities S and Q are both 
time dependent and, importantly, dependent on the solution so constructed. Only after these constraints are enforced 
would one actually have a solution. Pursuing this construction, although interesting, is outside the scope of the present 
paper and will be treated elsewhere. 

We observe that the equation ( l84l i is metriplectic. Indeed, by construction, we have a Poisson bracket (ITTI i (the 
Gardner bracket) and a symmetric bracket ( [83T l. Since these were constructed out of triple brackets, property (iii) of 
Definition in Section|4T|holds. Positive semidefiniteness of the symmetric bracket follows from dSTl i. 

The nature of the dissipation of (|84] i is of particular interest in that it involves the global quantities S and Q. This is 
reminiscent of collision operators, such as that due to Boltzmann a nd generalized no nlinear Fokker-Planck operators 



such as those due to Landau, Lenard-Balescu, and others (see, e.g., Morrison 1 1986ll ). The usual dissipation in 1 + 1 



systems is local in nature (see Sec. |4.5l l and dissipates energy. Thus the metriplectic construction of this section has 
pointed to a quite natural type of dynamical system that has dynamical versions of both the first and second laws of 
thermodynamics. The pathway for constructing other systems with nonlinear and dispersive Hamiltonian components, 
other kinds of dissipation, etc. is now cleared, and some will be considered in future publications. 



4.4.3 Some metriplectic generalizations 



It is evident that many generalizations are possible. We mention a few. 

• Without destroying the symmetries or formal metriplectic bracket properties we could allow one or both of the 
functions C and to depend on the field variable u or even contain pseudodifferential operations. In fact, such ideas 
were used in similar brackets in lFlierl and MorrisonI 1201 ill to facilitate numerical computation. 

• It is clear how to generalize (|64] l to preserve more constraints, say /i,/2,..., in addition to H. One simply first 
constructs the completely antisymmetric multilinear brackets {E,F,G,H,...} paralleling ( |63] l, and then, analogous 
to (|64|, constructs 



{F,G)haj2.-= / d0' / de"{Uie'),H,h,h, 
J51 Js' 



,F}^{e',e"){u{e"),H,iuh,...,G} 



(85) 



The bracket {F,G)hjiJ2,... guaranteed to be symmetric, conserve the invariants, and be positive semidefinite. 
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• It is of general interest to have metriplectic systems of the form 



F=^{H,F,G} + J^d9' J d9"{U{9'),S,F}'^{d',9"){Ui9"),S,G} 

(such as our example of Sec. l4.4.2] i for a suitably chosen function G; here H is the Hamiltonian and S is the entropy. 
Exploring the mathematics of when this is possible is an area to pursue. 
• The constructio n here is easily extendable to higher sp atial dimensions. For example consider the following triple 
bracket given in Bialynicki-Birula and MorrisonI 1 1991 ] : 



{E,F,G}^ f d''zEf[Ff,Gf] 

J Si 



(86) 



where z = {q,p) is a canonical six-dimensional phase space variable, f{z,t) is a phase space density, as in Vlasov 
theory, the 'inner' Poisson bracket is defined by 

[f,8]^fq-8p-fp-8c,- (87) 

We assume that the domain & with boundary conditions enables us to set all surface terms obtained by integrations 
by parts to zero, thereby assuring complete antisymmetry. Inserting the quadratic Casimir C2 ■— d^zf^ /2 into ( |86] | 
gives 

{F,G}vP = {C2,F,G}^ f d^zf[Ff,Gf] , 

the Lie-Poisson bracket for the Vlasov-Poisson system, as given in iMorrisonI il980ll . Thus, this bracket with the 
quadratic Casimir is formally akin to the construction given in Sec. 14.2.11 (although we note it reduces to a good 
bracket for any Casimir and in this way is like the case of so(3) of Sec. l4.2.2] i. The triple bracket of ( |86] | can be used in 
a generahzation of the bracket of ( l64l i to obtain a variety of energy conserving colhsion operators, with a wide choice 
of Casimirs as entropies. 



4.5 Hybrid dissipative structures 

Even if a system is not metriplectic, it is of interest to see if it can be obtained from an equation which consists of a 
Hamiltonian part and a gradient part with respect to a suitable Poisson bracket and metric, respectively. 

For KdV-hke equations, energy (the Hamiltonian) is generally not conserved when dissipation is added to the sys- 
tem. This is common for physical systems, but a more complete model would conserve energy while accounting for 



heat loss, i.e., entropy production. In the terminology of IMorrisonI 11200911 . models that lose energy, such as those treated 



here and those described by the double bracket formahsm of §27T| are incomplete, while those that do represent dynam- 
ical models of the laws of thermodynamics, such as metriplectic systems, are termed complete. Although incomplete 
systems do not conserve energy, they may conserve other invariants, and building this in, represents an advantage of 
various bracket formulations. Thus, we construct incomplete hybrid Hamiltonian and dissipative dynamics by combin- 
ing a Hamiltonian and a gradient vector field according to the prescription 

u, = {u,H} + {u,S) (88) 

where u n- {m,//} is a Hamiltonian vector field generated by H and u H> (m,^) is a gradient vector field generated by 
S (which could be H). Thus, ( , ) is, up to a sign, an inner product on the space of functions u. 

Consider the following examples: 

• With the usual KdV Hamiltonian of dSTT l and the Gardner bracket of (|30] | describing the Hamiltonian vector field, 
together with the choice 

1 P7L 

\2 



4-n J -71 

we obtain for the gradients of Corollary |2] 
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(i) M, = {u,H} — V^Hi = —iigee +6uug — u 

(ii) Ut = {u,H}-'VHi = -ueee+^uue + uee 

(iii) Ut = {u,H} - V^Hi = -ueee +(3UU0 ~ Jif{uQ) 

which is the KdV equation of ( [29l ) with the inclusion of a new term that describes dissip ation. Case (i) corresp onds to 
simple linear damping, case (ii) to 'viscous' diffusion, and case (iii) to the equation of Ott and Sudan 1 1969ll which 



adds a term to the KdV equation that describes Landau damping. For these systems the KdV invariant j^^Adu^ 
serves as a Lyapunov function. 

Choosing H = S = H\, the Kahler Hamiltonian flow of ( iMt together with the dissipative flow generated by ( 1211 1. 
yields 

Mr = {u,H\} — V^//i = — Mg — 

which describes simple advection with Landau damping. This equation possesses the damped traveling wave solu- 
tion. 

We note that we can derive the heat equation from a symmetric bracket of the form (|64] |, again with 5^(0', 9") = 
5(0' - 9"). Using this ^ and noting {U{9),Ho,F} = G'J9), we obtain 



iF,G)Ho^ idOF^G',. (89) 



Let us compute, for example, F{u) = {F, —H2)ho (see (l72b). Since 5H2/ 8u = — m, we obtain 

d9F,u^^F{u)^{F,H2)Ho^- [ d9F;,u'=[ d9 F.u" . 
at J51 J51 

This yields 

Ut = Uxx 

which is the heat equation. 

From these examples it is clear how a variety of hybrid Hamiltonian and dissipative flows can be constructed from the 
machinery we have developed. For example, if we replace the KdV Hamiltonian by H{u) = /^idS {^iij^{uq) + |m^) 
we obtain t he Benjamin-Ono equat i on wi th the various dissipative terms. Related ideas apply to fluid dynamics may 



be found in lGav-Balmaz and Holml 11201211 
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